Cosmology of hybrid metric-Palatini /(X)-gravity 
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A new class of modified theories of gravity, consisting of the superposition of the metric Einstein- 
Hilbert Lagrangian with an f(lZ) term constructed a la Palatini was proposed recently. The dynami- 
' cally equivalent scalar-tensor representation of the model was also formulated, and it was shown that 

, even if the scalar field is very light, the theory passes the Solar System observational constraints. 

Therefore the model predicts the existence of a long-range scalar field, modifying the cosmological 
and galactic dynamics. An explicit model that passes the local tests and leads to cosmic accelera- 
tion was also obtained. In the present work, it is shown that the theory can be also formulated in 
terms of the quantity X = k 2 T + R, where T and R are the traces of the stress-energy and Ricci 
£C) • tensors, respectively. The variable X represents the deviation with respect to the field equation 

trace of general relativity. The cosmological applications of this hybrid metric-Palatini gravitational 
theory are also explored, and cosmological solutions coming from the scalar-tensor representation 
of /(X)-gravity are presented. Criteria to obtain cosmic acceleration are discussed and the field 
equations are analyzed as a dynamical system. Several classes of dynamical cosmological solutions, 
depending on the functional form of the effective scalar field potential, describing both accelerat- 
ing and decelerating Universes are explicitly obtained. Furthermore, the cosmological perturbation 
equations are derived and applied to uncover the nature of the propagating scalar degree of freedom 
and the signatures these models predict in the large-scale structure. 
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I. INTRODUCTION 
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Recently, a novel approach to modified theories of gravity was presented, consisting of adding to the Einstein- Hilbert 
t-H ■ Lagrangian an f(7V) term constructed a la Palatini Using the respective dynamically equivalent scalar-tensor 
^ | representation, it was shown that the theory can pass the Solar System observational constraints even if the scalar 
field is very light. This implies the existence of a long-range scalar field, which is able to modify the cosmological 
and galactic dynamics, but leaves the Solar System unaffected. An explicit model that besides passing the local tests 
leads to cosmic acceleration was also presented. The main motivations to this approach lies on the fact that metric 
and geodesic structures of gravity could play important roles in dynamics as clarified by Ehlers-Pirani- Schild (see @ 
and reference therein for a comprehensive discussion on this problem). Furthermore, the good experimental results 
of general relativity, at local scales and in the weak field approximation, should be retained by any theory of gravity 
aimed to explain large scale and cosmological dynamics as prescribed by for example the chameleon mechanism |3| or 
by the disformal screening mechanism [4J. 

In this context, there has been a considerable interest in the modifications of the geometric part of the Einstein's 
field equations, motivated mainly by the observation of the late-time accelerated expansion of the Universe 
In particular, gravitational actions consisting of more general combinations of curvature invariants than the pure 
Einstein-Hilbert term have been investigated extensively [(| 0] ■ Einstein himself was more satisfied with the geometric 
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part of his equations, and has been quoted to say that while the left hand side is carved of marble, the right hand 
side is made by straw. However, in generalized gravity theories, the problem of coupling matter to gravity is often 
reduced to the question of in which frame matter resides with respect to gravity. The matter Lagrangian, and the 
corresponding stress-energy tensor, are defined in the usual way, but the metric that matter couples to can be related 
to the gravitational metric by a conformal, or, more generally, by a disformal transformation. This is a fundamental 
issue in relation to observations, essentially based on matter dynamics. Specifically, if one is adopting a Jordan frame 
matter is minimally coupled to geometry. On the other hand, in the so called Einstein frame, matter can result 
non-minimally coupled to geometry and the gravitational coupling can depend on time and distances Q. At least 
in the case of non-universal coupling, this fact can pose the fundamental problems of what the true physical frame 



non-conservation terms in the continuity equations, the general structure of the theory is retained. 

New light on this old problem is shed within the recently introduced framework denoted C-theories [l(| , that unifies 
the first and second order formalisms for gravity theories. In C-theories, the connection is related to the conformally 
scaled metric g^ v — C(7Z)g^ v with a scaling that depends on the scalar curvature 1Z. With nonlinear f(lZ), C-theories 
interpolate and extrapolate the Einstein and Palatini cases, and can avoid some of their conceptual and observational 
problems. It is important to stress that a bi-metric structure naturally comes out in this context by solving the 
Palatini equation for the connection [ll|, Ell an d the problem of the dark side of the Universe could be re-conducted 
to the confrontation of a metric in the Einstein frame to another one in the Jordan frame [l3| . 

In an earlier work [l4j . the known equivalence between higher order gravity theories and scalar tensor theories 
was generalized to a new class of theories. More specifically, in the context of the Palatini formalism, where the 
metric and connection are treated as independent variables (see [l5| for a recent review), the Lagrangian density was 
generalized to a function of the Ricci scalar, computed from the metric, and a second Ricci scalar, computed from 
the connection. These theories can be written as tensor-multi-scalar theories, with two or more scalar fields. It is 
important to stress that such an approach select the whole budget of degrees of freedom due to the gravitational field 
in generalized theories of gravity. Essentially, dynamics of any scalar field is described by an effective Klein- Gordon 
equation which comes out beside the standard Einstein field equations. In particular, a fourth-order metric theory 
of gravity can be described as a second order Einstein theory plus a second order Klein-Gordon theory; a sixth-order 
theory is equivalent to Einstein plus two Klein-Gordon equations and so on. The results is that effective scalar field 
description can put in evidence all the gravitational degrees of freedom of a given theory 1 Q. 

More radically, one may modify the response of matter to gravity by defining an action which depends nonlinearly 
upon the matter Lagrangian [rjj, or even its trace [l8j]. Generally, the motion is non-geodesic, and in fact, in these 
cases, the motion of matter is typically altered already in flat Minkowski space, and one may expect instabilities due 
to new nonlinear interactions within the matter sector. A natural way to obtain solely gravitational modifications 
of the behavior of matter emerges in the Palatini formulation of extended gravity actions. In this approach, the 
relation between the independent connection and the metric turns out to depend upon the trace of the matter stress 
energy tensor in such a way that the field equations effectively feature extra terms, given by the matter content. 
However, since the extra terms contain fourth order derivatives, the theory is problematical both at the theoretical 
and phenomenological levels [19|, [20| . 

In this paper, we will consider cosmological aspects of the hybrid metric-Palatini gravity considered in Ref. 
In particular, starting from the field equations derived assuming as variable X = k 2 T + R, comprehensive of matter 
stress-energy tensor trace and the Ricci scalar, we obtain cosmological solutions, where matter and curvature can be 
discussed under the same footing. 

The present paper is organized as follows. In Sec. UH we derive the field equations for f(X) gravity in D-dimensions 
showing the general aspects by which metric and Palatini features mix by using X as leading variable. Then we 
specialize to D = 4 and derive the cosmological equations. Standard matter and curvature scalar can be used to 
derive effective pressure and matter-energy density. The scalar-tensor representation of the theory is considered in 
Sec. Ill CI Here we show that /(X)-gravity can be recast as generalized Brans-Dicke theory but the fact that we 
are using a hybrid metric-Palatini theory has important physical consequences that allow to avoid shortcomings that 
emerge if one takes into account pure metric or Palatini approaches. In Sec. IHII we derive the cosmological background 
equations, solve them in an example model and formulate the model as a dynamical system. Section ITVl is devoted to 
the analysis of cosmological models emerging from /(X)-gravity in AD. We also go further to consider fluctuations 
around the background solution in Sec. |Vj where we derive the cosmological perturbation equations and use them to 
show that the structure formation is free of instabilities, with however some new signatures that can be used to probe 



Recent progress shows that a nonperturbative approach, where one takes into account an infinite number of higher derivatives, results 
in qualitatively different theory that may even exhibit asymptotic freedom for gravity [161 ]. 
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and test the models emerging from the new class of f(X) gravity theories. Conclusions are drawn in Sec. I VII 



II. /pO-GRAVITY 
A. The field equations in D-dimensions 



Let us start by considering the action for hybrid metric-Palatini gravity in /^-dimensions given by 



S = I d D s 



-g[R + f(K)+2K 2 £ m ] 



(1) 



where, in addition to the Einstein-Hilbert term and the matter Lagrangian which we assume to have the standard 
form, there is an extra term depending on both the metric and an independent connection T through 

7Z = g^lZ^u = g^ v {v " vol — f" ia u + f " A f * v — f ^x^av) ■ (2) 

The definition of the metric Ricci scalar R is, as usual, the above formula unhatted. By solving the equation of motion 
for the connection, one finds that it is compatible with the metric F(1Z) D ~ 2 g^ that is conformally related to g^, 
when the conformal factor is given by 

df(K) 



F(K) 



dlZ 



This implies that the Palatini Ricci tensor is 

D - 1 1 



F(TZ) <fl F(TZ) <v - -L-V^F(K)^ - ' 



D-2F 2 (K) v v F(K) 
Varying the action ((T|) with respect to the metric, we obtain 



(D - 2) F(K) 



g MV UF(H). 



+ F(K)K^ - \f{H)g^ = k 2 T^ , 

where the matter stress energy tensor is defined as usual, 

2 6(V=g-£ m ) 



We can solve 1Z from the trace of this field equation which yields 



D 



f(K) - F(K)K = -k 2 T- 



D 
~2 



1)R=X. 



(3) 
(4) 
(5) 
(6) 
(7) 



Assuming that the form of the function f(lZ) allows solution for the above equation, we can express 1Z algebraically 
in terms of X. The variable X thus measures the failure of the theory to satisfy the general relativistic trace equation 
which, for D — 4, gives R = —k 2 T. We can then express the field equation ([5]) in terms of the metric and X as 



Gfiv = -f{X)g tiy -F{X)R^ + F'(X)V tl X M 



1 



F"(X) 



D - 1 (F'(X)y 
D-2 F{X) 



D-2 

XnX M + K 2 T, lL , . 



F'(X)UX + F"(X){dXf 



(8) 



Note that (dX) 2 — X 4l X^. Due to the second order derivatives acting upon X, the theory is fourth order 
in derivatives of the metric. Note that here F(X) is the short-hand for F(7Z(X)), and F'(X) denotes F'(X) — 
dF(7Z(X))/dX, and so on. The trace of the field equations yields 



F'(x)nx 



F"{X) 



1 (F'(x)y 

2 F(X) 



{dXf 



D-2 



X + jf(X)-F(X)R 



0. 



Note that the Ricci scalar R, computed from the metric, is related to the Palatini 1Z via 

f i i / /* [ \ i \ ~ 

1Z{X) = R 



D 



(F'(X)Y 
\F(X)J 



OF(X) 



F(X) 



(9) 



(10) 



which is seen by contracting the gravitational field equation, given by Eq. @. Thus Eq. © and Eq. (fTU|) are 
redundant with Eq. ([7]). Let us now specify the theory in 4-dimensions, which will be explored throughout this work. 
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B. The gravitational field equations in 4-dimensions 



The above D-dimensional action can be specified as 

S = ^J d^xV^g [R + f(K)] + S„ 



(11) 



where S rn is the matter action, k 2 = 8irG, R is the Einstein-Hilbert term, 1Z = g^ v TZ^ v is the Palatini curvature, and 
TZ^v is defined in terms of an independent connection r"„ as 

TZ-nv = — r™a,r/ + r aA r A1[y — r MA r Q;y . (12) 

Varying the action given by Eq. (jlip with respect to the metric, one obtains the following gravitational field equation 



G\lv + F{K)K^ - -f(U)g^ = k% 



(13) 



The matter stress-energy tensor, as before, is defined by Eq. ([6]). Note that the independent connection is compatible 
with the metric F(TV)g^ v , conformal to g^, with the conformal 2 factor given by F(1Z) = df(7Z)/d7Z. This implies 
that 



= R^ + ^-E^rF^^Fin)^ 



1 -y^F(n). v -]--^—g^uF{n) 



2F 2 (TZ) y ''^ y F(K) v 2F(Ky 

The Palatini curvature, 1Z, can be obtained from the trace of the field equation (fT3"]) . which yields 

F(K)K - 2f(K) = k 2 T + R = X. 



(14) 



(15) 



As above, we can express 71 algebraically in terms of X if the form of f(7Z) allows analytic solutions. Again, the 
variable X measures how much the theory deviates from the general relativity trace equation R = —n 2 T. The field 
equation (fTB")) can be recast as G M „ = n 2 T^, where Tff = T^ v + T^ v , with T* v defined as 



1 1 



1% = ^l-f{X)g, w -F{X)R^ + F'{X)V^X. v + - F'(X)OX + F"(X) (dX)' 



2 F(X) 

The trace of the field equations is now 



X^X ;V 



F'(X)UX + 



F"{X) - 



1 (F'(X)Y 



(dXf + ±{X + 2f(X)-F(X)R}=0. 



2 F{X) 

while the relation between the metric scalar curvature R and the Palatini scalar curvature 7Z is 

3 



n{x) = r 



(F>(X)y DF(X) 
\F(X)J F(X) 



(16) 



(17) 



(18) 



which can be obtained by contracting Eq. (|14[) . Now we have all the ingredients to derive cosmological equations from 
this theory. 



2 If the gravitational lagrangian depended upon other curvature invariants besides TZ, one would obtain a disformal relation between the 
two metrics 
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C. Scalar-tensor representation of /(X)-gravity 

Like in the pure metric and Palatini cases [2l|, the action (|TT|) for f(X) theories can be turned into that of a 
scalar-tensor theory by introducing an auxiliary field A such that 

S = A / ^x^g [R + f(A) + f A (TZ - A)} + S m , (19) 



2k 2 j 

where /a = df/dA. Rearranging the terms and defining <f> = /a, V(<f>) — AJa — f(A), Eq. fp~9|) becomes 

1 
2~ 



S= — j d A xy^7j [R + <fiK- V{4>)\ + S m . (20) 



It is important to stress that in this case we are also considering a hybrid metric-Palatini theory. Variation of this 
action with respect to the metric, the scalar (f> and the connection leads to the field equations 

R^ + ^u- -(R+4>n-V) gflu = k 2 T^, (21) 

K-V+ = 0, (22) 
V a (V=ff&r) = 0, (23) 

respectively. 

The solution of Eq. ()23|) implies that the independent connection is the Levi-Civita connection of a metric h^ v = 
4>9tiv- This means that we are dealing with a bi- metric theory and IZ^v and R^ v are related by 

which can be used in the action (f20|) to get rid of the independent connection and obtain the following scalar-tensor 
representation that belongs to the "algebraic family of scalar-tensor theories" (22[, so that we finally arrive at the 
following action 



Hp, = R„v + ^d„4>d u( f> - - ( V M V^ + -g^dtf, ) , (24) 



1 

2^2 



d 4 x, 



-9 



(i + <£)R + ^ 9 m<W - V{<f>) 



+ S m . (25) 



This action differs from the w — —3/2 Brans-Dicke theory in the coupling of the scalar to the curvature, which in 
the w = —3/2 is of the form <pR. As we will see, this simple modification will have important physical consequences. 
Using Eq. (|2"4"|) and Eq. (|2"2"|) in Eq. (|21l) . the metric field equation can be written as 

(1 + </))iV = k 2 [T^ - + i 5/a , (V + n#) + V M V^ - ^d^d^ , (26) 

from which it follows that the spacetime curvature is generated by both the matter and the scalar field. The scalar 
field equation can be manipulated in two different ways that illustrate how this theory is related with the w = 
and w = —3/2 Brans-Dicke theories, which corresponds to the metric and Palatini scalar-tensor representations of 
/(-R)-gravity 0, respectively. Tracing Eq. ([21]) with , we find -R - $R + 2V = n 2 T, and using Eq. ([22]), takes 
the following form 

2V - (jjV^, = k 2 T + R . (27) 

Similarly as in the Palatini (w = —3/2) case, this equation tells us that the field <fi can be expressed as an algebraic 
function of the scalar X = k T + R, i.e., <p = 4>(X). In the pure Palatini case, however, <j) is just a function of T. 
The right-hand side of Eq. (|26p , therefore, besides containing new matter terms associated with the trace T and its 
derivatives, also contains the curvature R and its derivatives. Thus, this theory can be seen as a higher-derivative 
theory in both matter and metric fields. However, such an interpretation can be avoided if R is replaced in Eq. (|27[) 
with the relation R = TZ+ ^Ocf>— ^d^d^cj) together with TZ = V<p. One then finds that the scalar field is governed 
by the second-order evolution equation 



□0 + JLa ( ,08'>+— - - =?=-T, (28) 



which is an effective Klein-Gordon equation. This last expression shows that, unlike in the Palatini (w — —3/2) case, 
the scalar field is dynamical. The theory is therefore not affected by the microscopic instabilities that arise in Palatini 
models with infrared corrections [Tsl |20| . 

With these considerations in mind, we shall consider /(X)-cosmology in the scalar-tensor representation below. 
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III. /pO-COSMOLOGICAL FIELD EQUATIONS 
A. f(X) gravity and the late-time cosmic acceleration 



As a specific example of f(X) cosmological dynamics, let us consider the spatially flat Friedman- Robertson- Walker 
(FRW) metric given by the metric element 



ds 2 



-dt 2 + a 2 (t)dx 2 , 



(29) 



where a(t) is the scale factor. The Ricci scalar is given by R — 6(2H 2 + H), where H = a(t)/a(t) is the Hubble 
parameter, and the overdot denotes a derivative with respect to the cosmic time. As usual, the standard matter 
conservation law for the energy density p m and pressure p m is given by p m + 3H(p m + p m ) = 0. 
Thus, the modified Friedmann equations take the following form: 



K 2 p e s 



3H 2 



n 2 p cfi = -(2H + 3H 2 ), 

where p c g and p e ff are the total effective energy density and pressure, respectively, and are given by 



Pes = P j 

K 



>>-! 



F"{X) 



(F'(X)) 



21 



F(X) 



3 -F" 
2 



X 2 - ^F'(X) (X + XH) - 3F(X) (H + H 2 



)} 



Pes - p+^\\}{X)-^F'\X)X 2 -^F\X)(x + hXH) -F{X)(H + 3H 2 ) }, 



(30) 
(31) 



(32) 
(33) 



respectively. For an effective equation of state defined by the parameter ui e ff = Pes/peS-, we obtain w c ff = — 1 - 
2H/3H 2 . The late-time cosmic acceleration occurs when the strong energy condition (SEC) is violated, i.e., p e g - 
3p c ff < or u> c ff < —1/3. For simplicity, considering the specific case of vacuum p m = p m = 0, the SEC is given by 



PcS + <$PeS 



f(X)-6H F'{X)X + FH 



2 F(X) 



(34) 



As a toy model, let us consider the power law models f(7Z) = otR, n , where a and n constant parameters. Using 
Eq. (fl"5)) , we obtain 



f(X) 
F(X) 

F'(X) 

F"(X) 



X 



n-2 



X 



a(n - 2) 



n-1 
n-2 



X 



a(n - 2) 



-("-!) 
an(n — 2) 2 



X 



a(n - 2) 



As in the pure f(TZ) or "Palatini" gravity, the case n = 2 is degenerate, and in the following we assume n ^ 2. The 
general condition for the acceleration of the Universe can be formulated as 



1 



36H 2 n 2 X 2 + 9 (n - 1) (n - 1 + AHnX) X > -nX (6X + 2n - 4) 



X 



(n - 2) a 



l/n 



(35) 



To look for power-law solutions for the scale factor, a ~ t 13 , we plug such an ansatz into the Eqs. (|32 p -([33 ]l . 

For generic n, we find that the consistent solutions are described by the effective equation of state parameter, 
w c ff = Pos/PeS, given by 



Wcff 



Aj5 [1+ 4/3(1 - /3) - n/3(4/3 2 



+ 5)] 



f{n, a, 13, t) + Qnf3 [n/3(4/3 3 - 12/3 2 + 13/3 - 6) + /3(4^ 2 - 8/3 + 5) + n - 1] 



-1+3/5, (36) 
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where, for notational simplicity, f(n,a,/3,t) is defined as 

f(n, a, P, t) = na(n - 2)(2/3 - l)t 4 6 1/r 



/3(2/3 - 1) 



t 2 a(n - 2) 



(37) 



The first equality in Eq. (I36[) follows from inserting the ansatz in the field equations, the second one directly from 
the ansatz. One may now take specific values for the parameters a and n to solve away (3 and obtain u> c ff < — 1/3 to 
predict an accelerated behavior of the Hubble fluid. 

As a specific example, for the case of n = —1, so that the effective equation of state parameter, u> c ff simplifies to 



Wcff 



6/3 3 - 25/3 2 + 31/3 - 10 - | 



2 

-1+3/3. 



6/3 3 - 21^ 2 + 21(3 - 6 - ^t 4 
One of the solutions is the de Sitter expansion, 

Wcff = -1 , 

thus we expect that in the / ~ 1/1Z model the universe ends up in an accelerating phase. 



(38) 



(39) 



B. General dynamical system analysis 

To explore the general dynamical analysis, consider the /(X)-gravity action given by 

s = i / ^ x ^~3 ^ aR + fw\ + Sm 



(40) 



which is equivalent to the scalar tensor theory of the so called algebraic class 3 . Note that we include the parameter 
here for generality, and then the result in the previous section becomes 



s = i / d^V^a 



)R+^{d<pf-2K>v^) 



+ s„ 



where 



«M0) = I [r(0)^ - f{r{4>))] , r(<t>) = f \<f) 



(41) 



(42) 



The Friedmann equations can always be written in terms of the effective energy density and pressure, respectively 
defined as 



2 2 K 
3H = K PeS ~ 



(43) 
(44) 



k 1 K 

H = — (jOefif + Poff) + -o • 

2 a z 

For simplicity, we consider the flat universe K = in the following. For the theory provided by the action (j40|) we 
obtain the following modified Friedmann equations 



(Sl A + <t>)K 2 PeS = ~ — <f> 2 + K 2 V{(l))~3H^ + K 2 p mi 
40 

3 . ... 

(n A + 0) K 2 Pcff = - — <j) 2 -K 2 V{(f>)+(j) + 2H(j)+K 2 p r , 



(45) 
(46) 



respectively. 



3 By algebraic we mean that there always exists a functional relation </> = (f>(Q A R-\-K 2 T), but we stress that the scalar field is nevertheless 
dynamical except in the special case Q A = where one recovers the Palatini-/(7?,) theory |22l . 
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The conservation equations for the matter component and the scalar field are 

f)rn + 3H(p m + p m ) = , (47) 

4> + 3H^~^ + ^R- ~K 2 cj>V'(cj>) = . (48) 
Recalling that R — 6(2H 2 + H) and using Eqs. (l4"3")l and (|4"4")l . we can rewrite the Klein-Gordon equation as 

4> + 3H<p-^ + U'(<t>) + ^(p m - 3p m ) = , (49) 
where for notational simplicity, U'(cj>) is defined by 

u'{<t>) = [2V(4>) - (Oi + 0) v'm . (so) 

As a consistency check one can verify that the Klein-Gordon equation together with the matter conservation allows 
to derive flU) from (|35)l. By combining Eqs. 05} and (01]), we find that 

2V(<f>) - V\<t>)4> = ~ {VaR + K 2 T rn ) = hi . (51) 

The solution for <fi = <fr(X = 0) gives us the natural initial condition for the field in the early universe. The asymptotic 
value of the field in the far future may then be deduced by studying the minima of the function U (4>) defined by Eq. 

In order to study the dynamical system, we introduce the dimensionless variables 

r> — k2 Pm _ , k 2 V . . 

where N — log a is the e-folding time. The Friedmann equation ([4"5]) can then be rewritten as 

y 2 

Qa + x + y — z + — = Q m . (53) 
4x 

Due to this constraint, the number of independent degrees of freedom is three instead of four. We choose to span our 
phase space by the triplet {x, y, z}. The autonomous system of equations for them reads as 

%,N = V, (54) 
y.N = g^ AX { {3w m - 1) y 2 + Ax [{2>w m - 1) V - 3 (1 + w m ) z] 

-Ax 2 (1 - 3w m - 2u(x)z) + Att A [3x (w m - 1) y + y 2 - x 2 (2 - 6w m - Au(x)z)] | , (55) 

Z,N = ^ { ( 3w ™ _ 1) V 2 + 4a; [( 3w m - 1) y - 3 (1 + w m) z] 

+4f^£ (3 + 3w m + u(x)y) + Ax 2 (3w m - 1 + 2u(x)z) } , (56) 

respectively. We have defined u(x) = V'(4>)/V(4>). The relevant fixed points appear in this system. In particular, we 
have the matter dominated fixed point where x = y — z = and w e g = w m , and the de Sitter fixed point 4 that is 
described by w G ff = — 1 and 

x» = (2 — Hau^/u* , 2/*=0, z*=2/u*. (57) 

We denote the asymptotic values corresponding to this fixed point by a subscript star. In particular, the asymptotic 
value of the field x* is solved from the first equation in (|57|) once the form of the potential is given. As expected, this 
value corresponds to minimum of the effective potential (|50p. U'(x*) = 0. To construct a viable model, the potential 
should be such that we meet the two requirements: 



4 In addition, there exists the fixed point x = 



— corresponding to some kind of singular evolution. 
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• The matter dominated fixed point should be a saddle point, the de Sitter fixed point an attractor. Then we 
naturally obtain a transition to acceleration following standard cosmological evolution. 

• At the present epoch the field value should be sufficiently close to zero. Then we avoid conflict with the Solar 
system tests of gravity. 

Note that the simplest metric f(R) theories that provide acceleration utterly fail in both predicting a viable structure 
formation era and the Solar system as we observe it. The Palatini-/(7?.) models on the other hand can be ruled out 
as dark energy alternative by considering their structure formation or implications to microphysics, if such a theory 
is regarded consistent in the first place. As shown here and explored further below, f(X) gravity models exist that 
are free of these problems. 

To summarize: the field goes from </>j to 0*, where the former is given by 2V((f>i) = V'((j)i)(j>i and the latter 
by 2V(4>*) = (Ha + <t>*)V'(4>^). We just need a suitable function V((/>), i.e. f(lZ) in such a way that the slope will be 
downwards and (/>* near the origin. 



IV. COSMOLOGY OF /(X)-GRAVITY: ANALYTICAL SOLUTIONS 

As we have seen, the above theory of /(A)-gravity can be recast in terms of a non- minimally coupled scalar-tensor 
gravity. The scalar field represents the "deviation" of the theory from the standard general relativity where standard 
matter is also considered. In particular, the auxiliary variable X states "how much" the trace equation of the theory 
deviates from the general relativistic one. An important point to stress is that such a scalar field has a purely geometric 
origin and describes further degrees of freedom of the gravitational field coming from extended theories of gravity. 

It is straightforward to rewrite the cosmological equations in the absence of standard matter as 

H 2 = yPcff, (58) 

K 2 

H = -— (Pefif +Pcff), (59) 



where 



(! + (/>) K 2 p eS = -^-<j) 2 + K 2 V((j))-3H(j>, (60) 



(l + £)K>eff = --7T lf>-^V((f>) + 2H<P + (j>. (61) 
4</> 

It is important to stress that there is no standard matter in the definition of pressure and energy density, however, 
due to the definition of the scalar field cf> and its dependence on X curvature behaves as a perfect fluid. The scalar 
field satisfies the Klein-Gordon equation, 

4> + ZH4> - + 2 ~Y<t> W (</>) - (l + <p) V (0)] = o. (62) 

Note that as an indicator of the accelerated expansion one can consider the behavior of the deceleration parameter, 
given by 

d 1 H 3p off +p cff 

q= jtH- l = -ip- l = 2^r- 1 - (63) 

where accelerated expansion occurs when q < 0. 
From Eqs. (|60|) and (|6Tj) . we obtain 

{l + <f>)K 2 (p eS +p eS ) = -^-j> 2 -Hj> + 4>. (64) 

2(p 

By eliminating (j> with the help of the Klein-Gordon Eq. (f6"2"j). we get 

(1 + 4>) k 2 (peff + Pcff) = - 4ff0 - ^0 [2V (0) - (1 + 0) V {cj))] . (65) 

<p 3 
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Therefore the cosmological equations become 



3£P = 



2H 



1 



1 



1 



-— 4? -3H<j> + K 2 V((t>) 



-U 2 + AH^jy + ^0 [2V (0) - (1 + 0) V (0)] 



The explicit dependence on the scalar field for q is 

_ 3 f ij?l<\> + 4H4> + 2n 2 c\> [2V (0) - (1 + 0) V (</»)] /3 ' 
q ~ 2 | -30 2 /40 + k 2 ^ (0) - 3iJ0 

Clearly, the dynamics of cosmological models can be classified according to Eq. 
conditions on Eq. (|68l) assign, in principle, the functional form of f(X). 



(66) 
(67) 

1 . (68) 
Reversing the argument, 



A. Marginally accelerating models 

A particular class of models, which may be called "marginally accelerating", are those satisfying the condition 



5 = 0, which gives the equation 

2 /0 + AH4> + 2k 2 4> [2V (</)) - (1 + 4>) V (</>)] /3 

and then 

2 



-3^ 2 /A(f> + k 2 V {(j>) -3#0 



0, 



r- + 2Hcb+ ^n 2 4> [3V (0) - (1 + 0) V {dp)] = 0. 



(69) 



(70) 



The simplest marginally accelerating model can be obtained by assuming that the potential satisfies the condition 

3V(0)-(1 + 0)V'(0) = O, (71) 

which yields the following solution 

V(4>) =V (l + <P) 3 , (72) 
with Vq an arbitrary constant of integration. We also obtain immediately 



The Friedmann equation (|65[) becomes 



16. 



H = 



k 2 Vq (1 



from which one obtains the solutions 



4>(t) 



exp 






1 — exp 





(73) 



(74) 



(75) 



and 



a(t) = ao < exp 



Ak\J^ (t-to) 



1/4 



1 



(76) 



The importance of the marginally accelerating models is due to the fact that since the deceleration parameter is 
a monotonically decreasing function of time, models that reached the value q = starting from decelerating states 
with q > will end in an accelerating state with q < 0, as in fact can be seen from Eq. (|76p . In this way we have 
established the existence of at least one accelerating solution in the /(X)-gravity model. Furthermore, these types of 
solutions acquire a physical meaning being the possible junction between any dust-dominated era (q > 0) and recent 
accelerated expansion (q < 0). In some sense, they are the turning point between an epoch of structure formation 
and dark energy. 
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B. Accelerating models 

1. General considerations 

To tackle, in general, the problem of accelerating models, let us simplify the mathematical formalism by introducing 
the new variable x = In a, dx/dt = a/a = H. Then we have 

d d dx da ^ d 
dt dx dt dx a dx 



and 



d 2 ^ d ^ ^ d 2 dx H ^ + H 2 ^ (78) 
dt 2 dx dx 2 dt dx dx 2 



In these variables, the hrst modified Friedmann equation (1551) becomes 

1 + + 0,* + 4>%m. 



H 



2 



hi 



3 

On the other hand, the second Friedmann equation (|67l) takes the form 

H 1 f 1 2 2k 2 [2V (0) - (1 + 0) V (0)] 



(79) 



(80) 



H 2 2(1 + 0) U ' ' 3 H 2 
and the deceleration parameter is provided by 

1 /l^ 2| ^ , ^ [2V^W - (WW! (1 + + 0,, + 2 j40) | i 



2(1 + 0) ^ y(0) 

This result allows to formulate the condition for accelerating expansion as 

[2 (1 + 0) + 0^/0 + 40^] 0^ ( _20^ 

1/ 20(i + 0)(i + + 0, a + 0y40) 1 + 0' 1 j 

where we have used dV/dcf) — (dV/dx) (dx/dcf)) — (dV/dx) (l/0 iX ). The Klein-Gordon equation (|6"2")) can be trans- 
formed as: 

^ j.™ i 2^ 2 0[2T/(0)- (1 + 0)^(0)] 

0.XX + 30,2 - — i2 . + 0,^ + -g Jp ( 83 ' ) 

Now we search for explicit solutions that can be achieved by imposing the form of V(<j>). 

2. Power-law accelerating models 

In order to arrive at accelerating models, i.e., q — —H/H 2 — 1 < 0, we may assume that during the accelerating 
phase the deceleration parameter is a constant (in a pure de Sitter expansion the deceleration parameter is —1), so 
that 

q = — -rpr — 1 = —qo = constant . (84) 



This expression can be simplified as 



§2 = lo - 1 , (85) 



which provides the form of the scale factor a(t). Furthermore, for simplicity, we impose on the potential the condition 



27(0) -(1 + 0) V"(0) = O, 



(86) 



that gives the solution 

V(</>) = V (1 + cj>) 2 . 
With these considerations, Eq. (|83)l reduces to 

<l>,xx + (2 + q ) <j>,x - 7^<i>% = 0- 

By denoting dtfi/dx = u, we obtain d 2 4>/dx 2 = du/dx = (du/d<ft) (dcfy/dx) = udu/dcf), so that Eq. (1881) becomes 

du „ 1 o 

u _ + (2 + 50 ).--^=0 ) 

or equivalently 

1 . . 

which yields the following solution 

u{cj>) = C l j4>-2{q + 2)4>. 
Restoring the previous variables, one arrives at 

d(t> 



Ci^-2(g + 2 

1 



dx. 



In Z - C 2 



2(^ + 2) (^_2(g o + 2)0/Ci)' 
On the other hand, Eq. ([55)1 can be written as 

d 1 



and then immediately integrated 



H= a - 1 



a (l-<7o)t' 
to provide the solution 

a(t) = a t 1_<Zo . 

In this way we have obtained a power-law accelerating expansion depending on the value of go- 

C. The general case 

1. Abel equation 

In the general case, Eq. (1551) can be written as 

<t>, xx + 30,* - ^% + \yT^ { + 4(f> ' x + /(<W (! + ^ + ^ + ^L/ 4 } + 

f{4>)(i + <t> + 4>, x + ct>y±<t>) =0, 

where we have denoted £/(</>) = and / (0) = 2<jS [2 — (1 + 0) U{4>)]. Eq. (JSTJ) can be written as 

A ± 3 (fMO-\* 30(/ + 2) + /-2 ^ 2 4 + / ^ 

+ 2 (/ + m * + 4^(1 + 0) + 80(1 + 0) ^ + (! + «/ = 0- 
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Introducing as above the variable u, we obtain 



and by dividing with it 3 , we find 

1 du 3 1 30(/ + 2) + /-2 1 4 + / ,,,,^,1 . rin _. 

u z d<p 2 40(1 + 0) « 80(1 + 0) u 01 

By denoting u = 1/v we obtain an Abel equation of the form 
Our task is now to integrate such an equation. 



2. A simple accelerating solution 
The simplest case corresponds to / = 0, V (0) — Vq (1 + 0) 2 , and <C 1. Then the Abel equation takes the form 



dv 1 o 1 / 

d^" 3 " -20=°' (1 ° 2) 



v = MVgt^. (103) 



with the general solution 



For small values of 0, the solution can be approximated as 

Since <C 1, one can approximate the potential as V ~ Vq = constant. Then the condition for accelerating 
expansion takes the form 

[2 + 0y0 + 40j0^ ^ 

20(i + 0, + 0y40) +2 ^ >0 ' (105) 

which is obviously satisfied since -x > if C > 0. 



D. Parameterizing by g 

In view of a cosmographic analysis of the models, it is convenient to parameterize the cosmological equations by 
the deceleration parameter q. As discussed in [23j j . cosmography is an extremely useful tool to discriminate among 
concurring cosmological models and, such a feature, is strictly related to suitable parameterizations. In our case, the 
approach is very straightforward, starting from Eqs. and (155)) . 

Thus, the field equations given by Eqs. (|66l) and (p7|) . in terms of the deceleration parameter, can be recast in the 
following extremely simple form 

3H 2 = ^^[2F(0)-0^(0)], (106) 

h = -'j^TW [2V (<w " w> m ' (107) 

where the deceleration parameter is defined according to the sign convention of Eq. (|84|). Therefore, the cosmological 
field equations can be expressed in terms of the scalar field, the scalar field potential, and the deceleration parameter. 
However, the Klein-Gordon equation, given by Eq. (|62p . also needs to be taken into account. Let us work out some 
simple examples adopting this parameterization. 
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1. Exponential expansion 

In the case of the exponential expansion q = — 1, H = Hq = constant, H — 0, and a(t) = ao exp [Ho (t — to)], 
respectively. For this case, Eq. (|107|) is automatically satisfied, while Eq. (|106[) becomes 



which provides the following solution 



cj>V' (<f>)-2V( ( f>) + -H 2 = 0, 



v{^) = \hI + ^, 



(108) 



(109) 



where Vo is an arbitrary constant of integration. The general solution for the scalar field is obtained from the 
Klein-Gordon equation, which is given by 



'(*) 



3!Zo 



exp [Ho (t - to)] ± V3Vo 



(110) 



The physical solution corresponds to the choice of the positive sign, in order to avoid any singular behavior of the 
scalar field. 



2. Power law accelerated expansion 
As above, in the case of power law accelerated expansion, we have the constant q = qo and 



H(t) = 



1 



H = —- 



1 1 



with qo =/= 1. 



The field equations, given by Eqs. (|106j) and (|107j) . provide the following expression 

(l-qo)^ 3(1 + go) 
Hence, the field equations reduce to a single condition for the potential, given by 

3 (1 + qo) 1 



2V (cp) - 4>V {4>) 



(l-9o) 2 t2 ' 



(111) 



(112) 



(113) 



Reversing the argument, any form of scalar field potential (and then the /(X)-function) compatible with power law 
expansion a(t) = aot 1_9 ° can be achieved in this way. 



V. COSMOLOGICAL PERTURBATIONS 



A. Field equations and conservation laws 



To understand the implications of these models in the context of structure formation, we will derive the perturbation 
equations and analyze them in some specific cases of interest. This paves the way for a detailed comparison of the 
predictions with the cosmological data on large scale structure and the cosmic microwave background. For generality, 
we will keep the parameter Qa in the formulas in this section. 

We work in the Newtonian gauge [24j |. which can be parameterized by the two gravitational potentials $ and ^, 



The — part of the field equations is 



(1 + 2*) dt 2 + a 2 (t) (1 + 2$) dx 2 



(114) 



k 2 
a 1 



3 \H 



2(0, 



$-3 F 



H<f) 



1 

2n A + . 



tt A + c/) 40 (tt A 
3 {H 



fc 2s 



a 



20 



(115) 
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where we have denoted tp — 8(f). The Raychaudhuri equation for the perturbations reads 



(h 2 



+ 2H 



6 



k 2 



a 2 Q A + 1 



H(j> 



1 



* + 3 \2H 



k 2 



2<M 



6$ 



The — i equation is 



k 2 (S Pm + 38 Pm ) + [6H 2 + 6H + 3^ - 2V'{4>) + —\ v> + 3[H--^\(p + 3ip 



H 



$ + $ 



i 



2 (ft A + 0) 



K 2 (Pm + Pm) av m + I H + ^ j ip + if 



Note that the set of perturbed field equations is completed by the off-diagonal spatial piece: 

* + $ = — . 

n A + <j> 

Assuming a perfect fluid, the continuity and Eulcr equations for the matter component are 



5 + 3H (c 2 - w) 5 = - (1 + w) (s® - ^-v^J , 



(116) 



(117) 



(118) 



(119) 



(120) 



respectively. The linear part of the Klein-Gordon equation is then compatible with the above system. For complete- 
ness, it is 



24> + 6H<j)- ^<j) 2 ^j * + - 3$) - ^SR. 



This completes the presentation of the field equations and the conservation laws. 



B. Matter dominated cosmology 

Let us consider formation of structure in the matter-dominated universe, where w m — c 2 — 0. Note that at the 
relevant scales, spatial gradients are more important than the time derivatives and, consequently, the matter density 
perturbations are much stronger than the gravitational potentials. Combining the continuity and the Euler equation 
at this limit, one obtains 



k 2 

-2HS-^. 

0T 



(121) 



We need then to solve the gravitational potential. Let us define II = a 2 p m S m /k 2 and write the field equations and 
the Klein-Gordon equation at this limit in a very simple way as 



(fi.4 + 0)$ = U-if, 

(n A + <t>)(y + $) = -tp, 

-2 (Q A + </>)* = U + ip, 

3<p = -20 (* + 2$) . 



(122) 
(123) 
(124) 
(125) 



We immediately see that one of the equations is (as expected) redundant, and that the potential * is (as usual) 
proportional to II, where now the proportionality is given as a function of the field (p. Our result is 



5 + 2H5 = 4irG cSPm 5 , 



with 



G. 



off 



n A -\cj> _ G 



(126) 



(127) 



n A (n A + (/)) 

This shows that there are no instabilities in the evolution of the matter inhomogeneities, in contrast to the Palatini- 
f(lZ) models and some matter-coupled scalar field models (recall our theory can be mapped into such in the Einstein 
frame) . 
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C. Vacuum fluctuations 

The propagation of our scalar degree of freedom in vacuum is also a crucial consistency check on the theory. Let 
us set p m = 0. Let us consider the curvature perturbation in the uniform-field gauge £. In terms of the Newtonian 
gauge perturbations this is 

C = ( 128 ) 

<P 

After considerably more tedious algebra than in the previous case, we obtain the exact (linear) evolution equation 

4> + 2H{n A + ct>) 1 2 \<t>(n A + <f>) i - ^ (n A + (j>) 3 J 

The friction term depends on the perturbation variable we consider, but the sound speed term retains its form. Thus 
perturbations at small scales propagate with the speed of light, as in canonic scalar field theory. This excludes also 
instabilities in the scalar field perturbations. Now equation (j!29[) can be used to study generation of fluctuations in 
/(X)-inflation. Construction of specific models and their observational tests are left for forthcoming studies. 

VI. CONCLUSIONS 

In this paper, we have discussed cosmological applications of hybrid metric-Palatini /(A)-gravity. This theory differs 
from Palatini /(72.)-gravity for two main aspects: i) both metric and Palatini approach are considered including R 
and f(TZ) in the gravitational action; u) deviations from standard general relativity are considered including the field 
X = k 2 T + R where X is exactly zero in the Einstein case. This approach is particularly relevant in view of the 
so-called chameleon mechanism useful to parameterize the astrophysical scales according to the cosmic densities: At 
short scales (Solar System scales) X — > 0, while it increases at larger scales since dark matter and dark energy effects 
have to be included. Another option is to consider the field dynamically evolving to X — > today, which renders 
cosmological dynamics compatible without resorting to the chameleon mechanism. Considering /(A)-gravity means 
not distinguishing a priori matter and geometric sources, but taking into account both of them into the dynamics 
under the same standard. Furthermore, shortcomings connected to the conformal transformations are avoided. For 
example, the singular cases w = for metric /(i?)-gravity and w — —3/2 for Palatini /(7£)-gravity in the Brans- 
Dicke-like representation of these theories are avoided, f{X) gravity being naturally metric-affine. This means that 
our hybrid model can be fully recast into a scalar-tensor theory of gravity, where the kinetic and potential terms are 
well defined. Furthermore, it is straightforward to recover general relativity, which is fully restored as soon as 4> — > 
(see Eq. (J25J). 

Here, we have also discussed the FRW cosmology coming out from /(A)-gravity. Accelerating expansion, and, in 
general, any cosmological behavior, strictly depends on the effective scalar field potential, which assigns the form of 
f(X). In particular, we have examined toy models where the deceleration parameter evolves according to such a 
potential. This approach is particularly useful in view of a cosmographic analysis to discriminate realistic models by 
observations (see [25| for /(i?)-gravity in metric approach). However, we conclude in noting that this is only a seminal 
paper. An extensive study of /(X)-cosmology and its comparison with observations, amongst other issues, will be 
the topic of forthcoming papers. 
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